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Abstract 

A model for three dimensional detonation is proposed based on the approxima- 
tion that the detonation thickness is small compared to the characteristic scales 
of the fluid motion. In this framework detonations are treated as a modified 
hydrodynamic discontinuity. The altered Rankine-Hugoniot jump conditions 
take into account the internal structure of the detonation including the chem- 
ical reaction. The position of the discontinuity surface and the corresponding 
jump conditions are derived from first principles. The final modified conditions 
are dependent on curvature, flame thickness and stretching and allow for simple 
physical interpretation. 
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Scale flow motion 
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1. Introduction 

Accidents involving the release and subsequent combustion of burnable gases 
can be classified, utilizing the flame propagation regime, ranging from slow de- 
flagrations to detonations. The latter appears to be particularly important since 
detonations are usually considered as a worst case, often resulting in accident 
scenarios with devastating consequences. 

During the last decades numerous modeling studies (e.g., pQ, [2], [3], [4]) 
showed that the detonation process in combustible gaseous mixtures can be 
successfully reproduced if the internal structure of the detonation wave is re- 
solved. Good practice implies resolving the 3D structure of detonation cells, 
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but ignoring details of the internal structure of the shock or of the chemistry 

H, ffl- 

The thickness of detonations is typically small compared to the character- 
istic scales of the fluid flow. Jin and Stewart [5] and Klein et al. 6| utilized 
the limit of asymptotically thin detonations to derive the asymptotic theory 
for weakly curved detonations. Thus, a simplified model for detonations can 
be built where the fuel consumption zone shrinks to an infinitely thin surface 
of discontinuity separating the cold mixture from the hot products. For de- 
flagrations, e.g. laminar premixed combustions, models considering flames as a 
gas-dynamic discontinuity exist since the pioneering works of Darrieus and Lan- 
dau who derived the jump conditions across the flame. More recently, Matalon 
and Matkowsky [7 considered arbitrary flame shapes for nearly equi-diffusional 
flames with thermal expansion in general flow fields. The leading terms of the 
jump conditions were those of the Darrieus-Landau model and perturbative cor- 
rections were obtained in the next order of approximation. Their formulation 
is based on high activation energy asymptotics. Klimcnko and Class [8] em- 
ployed tensor calculus and orthonormal coordinates to simplify the derivation 
of the flame speed relation of Matalon and Matkowsky. Their approach ex- 
plicitly exploits the distinctiveness of length scales instead of high activation 
energy, i.e. large Peclet number Pe, defined as the ratio of flow length scale 
to flame thickness. In the subsequent series of papers, their methodology was 
applied to multi-step chemistry [§] and wider reaction zones [TO] . Finally, Class 
et al. reconsidered the Rankine-Hugoniot jump conditions for the flow 

field and showed that perturbative correction of the jump conditions depend on 
perturbative corrections of the location where these jump conditions are evalu- 
ated. Moreover, they shown that there exists a unique location within the flame 
structure, where no extra inertia terms arise due to the continuity surface. This 
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corresponds to a vanishing perturbative correction of surface mass. 

In the present work, the authors apply the Class et al. [TT] methodology to 
detonation. Modified Rankine-Hugoniot jump conditions are derived, implicitly 
including the full effect of the chemistry. The methodology developed by Class 
ct al. [TT] has the essential advantage that perturbative corrections of the jump 
conditions are expressed as functions of the leading order solutions, so that the 
explicit evaluation of perturbative corrections of the solution becomes outdated. 
For detonation, this implies that the results of the Zeldovich-von Newmann- 
Doring theory can be utilized as a leading order planar model to close the 
system [JJ, |12j and obtain the final modified jump conditions which include 
perturbative correction terms. 



2. Analysis 

2.1. Governing equations 

The reactive Euler equations are considered, 



d t <t> + V-J(<l>) = Q(<t>), 



(1) 



where for the continuity, momentum, energy and species 
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(2) 



The system can be made dimensionless with the use of the reference variables p s , 
c si If where the sub index s represents the conditions at the Von Neumann peak 
and If is the characteristic scale of the flow motion. The dependent variables 
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can be written as p — pp Sl v — vc s , p — p ■ p s c^./j, e — ec 2 s , Y = XY S while the 
independent variables become t = t ■ If jc s , x = xlf. The Arrhenius chemical 
consumption rate can be written using the variables 9 a — E a /(R g T s ) and k = 
(kl c )/c s ) as 

W = k(l-Y)cxp(-0 a p/p). (3) 

We assume a large ratio Z of the hydrodynamic typical length to the half 
detonation thickness Z — If jl c . A thin detonation structure corresponds to 
an intense source term Q which is re-scaled accordingly Q — Z/Z ■ Q = ZQ' 
yielding 

d t <t> + V ■J{4>) = ZQ I {4>). (4) 

It is convenient to transform the equations to a moving generalized curvi- 
linear coordinate system. These coordinates were utilized in the references [S], 
01 [10], [II]- The coordinate system is attached to the discontinuity surface 
with its normal direction pointing towards the products. Its tangential direction 
moves with the local tangential flow. In this system, the flame is at rest with 
no flow along the flame surface. The tensor calculus utilized in our analysis is 
mainly contained in reference [13] . The coordinates [11] are orthonormal with 
x 1 direction perpendicular to the surface of the flame, so that x 1 coordinate 
lines are normal to the surfaces x 1 — const. The contra-variant metric tensor 
is given by g 1 ^ . Due to the orthogonality g la = 0, and due to normalization 
g 11 = 1. In the curvilinear coordinates, the system of differential equations can 
be written as 

d t (v50) + d xi (VgJ j (<!>)) - VgZQ'W (5) 
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with the flux vectors as J- 7 ((f)) = (V —u J )(f> with u J representing the speed of the 
moving coordinates relative to fixed Eulerian coordinates. The fluxes become 
according to [IT] . 
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. Decomposing the Eq 



(6) 



and tangential directions and introducing the stretched normal spatial variable, 
X = Zx 1 , Eq. (5) yields, 



z- 1 [d t (V9<t>) + d x . (VgJ a W)} + d x (V~9 'J 1 (*)) = VgQ'tt)- 



(7) 



2.2. Asymptotic derivation of the fluid equations and jump conditions 

Since Z, the ratio between the length of the fluid flow and the consump- 
tion zone thickness, is assumed to be asymptotically large, the variables are 
expressed in terms of an asymptotic series expansion in powers of l/Z , i.e. 



(0) 



Z~ 



(i) 



0(Z 1 ). The volume element is Taylor 



expanded around the discontinuity Jg = Jg 



(0) 



g. Z- 1 X + o(X 2 



Addi- 



tionally, we make use of the equalities y/dm — —2H and \ = dt(y/g ^) / y/g ^ 
with mean curvature H and stretch \- Equation ^ becomes 



z- 1 

+ d x 



(d t + x) (0(o) + ^V(d) + (g m r 1/2 d x * (v5 (0) (^o) (4>) + z-'Jw); 
1 - z- x mx) (4)^) + z- 1 4,)] = (t - z-Hex) (Q[ 0) (<t>) + z- 1 q' {1] ). 

(8) 
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Collecting coefficients of like powers in Z n , the zeroth, Z°, and first order Z 1 
terms are 



(dt + X) 0(0) + (5(0))- 1/2 ^ (Vfl(o)^o) W) + 

+ 9x (J ( \) - 2#XJ ( 1 0) (</,)) - Q' (1) (0) - 2ffXQ' (0) (</,). 



(9) 



(10) 



The leading order normal fluxes are, 
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(11) 



The fluxes in the momentum equation can be decomposed into the terms em- 
anating from the tangential and normal components. Using a parallel vector 
field li normal to the flame at the point of consideration, 



J( )(pv % h) = "i( 0) U( 0) Zi (0) +p ( o)/7 ' ' 



(o) ' 



(12) 
(13) 



The perturbative correction (first order) terms require some manipulation. 
The first order normal momentum flux is, 

J(i)(pv l k) = ( m (o) u (i) + m (i) v (0) +P(i)/7-5(o) +P(0)/7-ff(i)) + 

( m \o) v (o) + P(0)/7 ' ff(o)) (I 4 ) 



+ 



with fl n = 1, <$) = 1, 



9(1) =0 Vi^ and (mj 0) vf 0) + f> (0 )5$)) = 
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yielding, 



J(i)(pv%) = ( m ( ) u (i) + m (i) w (o) +P(i)/7'.9(o)) ii (0) - 



(15) 



JhJpe) simplifies substantially due to g %1 = 0, V i ^ 1. Ultimately, the first 
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(16) 



The asymptotic expansion should be also performed in the source term Q ss 
Q(0) + • The calculation involves asymptotic expansion of the variables 

inside and outside the exponential function plus a Taylor expansion of the latter. 
Finally, the first and second order terms are 



Q (0) = P (0 )fc(l - A (0) )e(-M p <°> T <°>)), (17) 
Q(D = (P(i)*(l - A (0) ) - Pi o)k\ (1) ) e (- fl »/(fco)T C o))) + 

+ [ P(0) fc(l - A (0) ) (T (1) /T (0) +P ( i)Mo)) e„/ (p(o)T ( o))] e^M^co,)). 

(18) 

#.3. Planar detonation 

The Eq. ([9]) combined with (11) and (17) constitute an equation system for 
planar detonations. This system is equal to the known equations utilized to 
derive the classical results of the ZND theory of detonation, see [12], [I], [14] . 
The ZND theory provides an analytic expression for the pressure, velocity and 
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specific volume profiles [H] as a function of the reaction progress variable A, 



p = a+(l-a)(l-b/3X)*, v = (1 - p)( 7 M s )- 1 + M s , T = v/M 8 . (19) 



The auxiliary variables appearing in Eq.(19l are, 



D = D/c s , M s = ^l^lly f> = Ql/cl 

2 7j D 2 -( 7 -1)' (l-a 2 )( 7 +l) V ' 



Equations (19) allow expressing the half reaction zone length as, 

l c = tsk- 1 !} v(X)(l - A)- W^ A > T ( A »dA. (21) 

It is interesting to note that the whole system can be converted to the spatial 
formulation using the change of variables, 



X/Z = x = f Q \(()/W(Od(, (22) 
dx/dX = v(X)/W(x)- (23) 



2.4. Modeling 

A detonation can be considered as a small-thickness layer separating the 
fresh mixtures of the burned products. We propose the derivation of a three 
dimensional hydrodynamic model in which the internal structure of the detona- 
tion as well as the chemical reaction is substituted by modified jump conditions. 
Conceptually, this construction asymptotically extends the simplest planar sta- 
tionary theory to three-dimensional non-stationary flow. 

In the derivation we consider two models for the detonation simultaneously, 
the hydrodynamic and the reactive model. Away from the consumption area, 
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these models are identical. The jump conditions of the hydrodynamic model, 
particularly the position and the amplitude of the discontinuity, are going to 
be determined from the internal structure of the detonation (reactive model). 
Therefore, two sets of equations must be handled each corresponding to one of 
the models. Capital letters are used to designate the reactive and lower-case 
for the hydrodynamic model. The system of equations ([9| and ( 10 ) may be the 
re-written as, 








Q(o)(*). 



(24) 







(dt + x) 



2HJ (o) 



0(0) , ,-i/ 2 q r ja ^(o) 
- (,9o ) d x « | v£ (0) J (0) 1 



(o). 



(0) 



(25) 



In the detonation profile, see Figure [T] both models coincide in the ini- 
tial, rarefaction and final states. The models exclusively differ in a thin zone 
surrounding the shock, i.e. between the hydrodynamic discontinuity and the 
Chapman- Jouguet point, see Figure [2] In the area of appreciable chemical re- 
action the hydrodynamic model is an extrapolation of the rarefaction wave. For 
simplicity, we set the origin of coordinates in the discontinuity of the hydro- 



dynamic model. Therefore, the Eq. (24 1 and (251 differ with respect to jump 



position, jump conditions and reaction source term. 



Figure 1. 



Figure 2. 
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Subtracting ( 24 ) from ( 25 ) and rearranging we obtain, 



4)W) = Q(o)> (26) 
9x(J ( i)(<I>) - J(\)(0)) = g' (1) - 2HXQ{ 0) (cj ) ) + d x (2HX (j^) - J ( 1 0) (<^))) 
- (9* + X) (<f ( o) - 0(0)) - {%T 1/2 d x « (yg (0) (j ( Q 0) ($) - J ( " o) (0))) • (27) 



Equations (26) and (27) can be piecewise integrated from — oo to oo. Taking 
into account the discontinuities and applying the Fundamental Calculus Theo- 
rem, 
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Jk)(4>) „ = (Q(i)(*) - 2ifXQ (0) ($)) dX- 



CJ 



+ a* (.r ; 2hx (ji 0) ($) - 4 } (0)) dx) - (ae + x) rr (*«» - * ( t») 
- rr (5o)^ l/2 ^ (v9 (0) (4o)( $ ) - J (o)(0 



(28) 



dX- 



(29) 



We have designated with the index VN the reactive discontinuity and with 
the index CJ the hydrodynamic, by analogy of its conditions with the CJ point. 
The discontinuity in the reactive model is a shock. Therefore, it is infinitely thin 
and the conditions applying for the shock discontinuity, see i.e. |15j , should also 
govern the curved case. Furthermore, the composition does not change across 
the discontinuity since the reaction starts at the high pressure side of the shock. 
Therefore, 



{[j 1 (*)] vn = o, [4>(*)] w = °}=> [ J w ($) 



= 0. 



VN 



(30) 



Substituting the condition JiL ($) = J^ 0) (cj)) for X — > — oo into Eq. (28) we 
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obtain for the vector 8 = (p, pv*li, pe) 



(31) 



This has several implications. The mass flux satisfies Ml 



(») 



(0) 



0. In 



the moving curvilinear coordinates, MfiL — Tn? \ — 0. Furthermore, the tan- 
gential component of momentum is continuous across the shock. Therefore, 
Jfa(RV%) = P( )/7 ' 9(o)k (0) and thus, defining II = (P (0) -p( ))/7 we obtain 
J" ^(RV l li) — J" - j (pv l li) = Hg l ^li {0) . Finally, the normal component of the mo- 
mentum flux, Jhy.(RV l li)l\ is identical in the reactive and CJ hydrodynamic 
model. Therefore, II = M^V^ — wjUwjL. 



The Eq. (29) can be simplified applying the previous considerations. For 



continuity, momentum and energy the first two terms in the RHS of ( 29 1 vanish 



Moreover, the last term of the RHS cancels for continuity due to vanishing 



transverse momentum. For the continuity Eq. (29) reduces to 



CJ 



(^ + x)L OO oo (i?[ )-P(0))^- 



(32) 



The simplifications in the momentum equation are, 



. m (0) v (l) + m \l) v \o) + 3(0)P(1)7 ) l Ho) 



Ho & + x) (fl(o)V? 0) - p (0 )«{ 0) )i i(o) dX, (33) 



those of the energy equation are, 



m (i) e (o) +W(o) e (i) +P(i)/l ■ 9(1) 'So)Mo) +P(o)h-9(o) ( m Vp) ( i) 
= -J-oo ( 9 t + X) ( R (o) E (o) - P(o)e (0 )) dX, 



CJ 



(34) 



13 



and finally those of the species equation yield, 



™(D A (°) + mJ 0) A (1) cj = (Q W (RL) - 2HXQ (0) (RL)) dX- 



(dt + X) (R(o) L (o) - P(o)A( )) dX. 



(35) 



We may define I R = (i? (0) - p (0) ) dX,I a = f™ x UdX and J E = ( E (Q) £ (0 ) - 
P(o)e(o))^- The equations, i2(o)V ( o) -P(o)«{ ) = "(ojCfyo) _ P(o)) and (3o) _1/2 ■ 
^"(^(0)^(0)^(0)) = ~ 2iJ/ l + (9oy 1/2 dx<*{y/g {0) g(o)h m )> see E3, allow ma- 



nipulating ( 32 ) , ( 33 ) and (34 1 to obtain 



'(i) 



CJ 



= -{d t + x)iR, 



(36) 



CJ 



( m lo) v (i) + m (i) u (o) +P{i)h - 5(0)) ii (0 
= - (fit + X ) w l (0) Z i(0 )/ii - (2^i + g af3 l p d xa ) I (j , (37) 

TO (i) e (o) + ™(o) e (i) +P(i)/7-ff(o) 'W(o)Mo) +P(o)/7'5(o) ( m Vp)(i)J CJ = 
= -(9t+x)/ E . (38) 



We combine the leading order jump conditions with the perturbative cor- 
rections and decompose the jump condition for the momentum in normal and 
tangential components. The result obtained is 
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Z^+oiZ- 1 ). 



(39) 



The first row of the Eq. ( 39 ) shows that the normal mass flux m 1 experiences 
0(Z~ 1 ) variations in the detonation structure. This represents the excess of 
mass in the hydrodynamic model compared to the reactive model. Following the 
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methodology described in [TT] , we define without loss of generality the position of 
the artificial discontinuity in the hydrodynamic model requiring identical normal 
mass flux in the fresh and burned mixtures. A continuous mass flux across the 
discontinuity surface is obtained. The condition [^ 1 ] CJ = requires Ir = 
where, 

Ir =/_«,(-%)) - P(o))dX + f* VN (R {0) - p {0) )dX+ 

+ 0^(0) - P(o))dX + C cj (R (0) - P (0) )dX. (40) 

The areas under the curves in Figure [T] and [2] delimited by the hydrodynamic 
discontinuity and the Chapman- Jouguet point cancel. Therefore, the first and 



last term of the RHS of Eq. (40) cancel. The difference of slopes, see Fig- 
ure [2j between rarefaction and consumption curves is asymptotically large. In 
stretched coordinates, the slope of the hydrodynamic model is zero at leading 
order and does not affect the integral Ir in the consumption area. After some 
manipulation we obtain, 

7 « ~ Ixvn ( R (P) ~ Pcj) dx ~ (pej - Po)X V n, (41) 

as can be graphically confirmed in the Figures [l] and [2j The position of the 
detonation shock relative to the artificial discontinuity is 

X V n « (PC J - A,)" 1 f£™ (i? (0 ) - Pcj) dX. (42) 
By similar considerations we find, 

^ = 7~7x (^(0) - P(0)) dX « 7- 1 (h ~ (pcj - Po) X VN ) , (43) 



where I\ = J^° J (P(o) ~ Pcj) dX. The terms inside integral I\ may be refor- 
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mulated changing the variables according to Eq.( 23 ). Now, the more convenient 



composition formulation is recovered, and thus I\ = Z (P(A) — Pcj) U (A) 1 



(r(A)) dX = ZI' X . Similarly Eq. (42 1 can be expressed in the composition 
formulation yielding, 



x VN 



Pcj ~ Po 



Pcj ~ Po 



Finally I a is reformulated, 



(44) 



I* =7 - (pcj - Po){pcj - Po) 1 I P )Z, 



(45) 



and Is is rescaled Is = I^Z. The Eq. (39) is rewritten as 



rr^v 1 + p/7 
m 1 vP 
m 1 e + m l p/ pj 



id t +x) 



cj 



M 







2HJ- 1 

g af) T l d x c 





Pcj - Po 
Pcj - Po 



(46) 



The integrals I\ and I p can easily be understood considering the Figure [2j 
Both integrals represent areas contained between the detonation curve and the 
horizontal CJ conditions. I\ and I p depend on the chemistry model. Even 
numerical evaluation for complex chemistry is possible. In the present work, the 



explicit expressions, Eq. (19), were selected in the analysis. Obviously, chemical 



conversion will not fully complete in a layer of finite thickness. This does not 
contradict our assumption of a consumption layer of small finite thickness, since 
the chemical conversion becomes exponentially small for A — >• 1. Accordingly 
the integrals I\ and I p converge for A — > 1. Note that, analogous to references 
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[TT] and [TB] the stretch x can be calculated as 



<o))\x=q+ 



(47) 



Here the stretch experiences a jump at the hydrodynamic discontinuity surface. 
Thus, the evaluation is taken on the high pressure side of the interface, a fact 



that is emphasized by the symbol + in Eq. ( 47 ) 



2. 5. Final jump conditions 

A clearer understanding of the jump conditions may be obtained re-writing 



Eq. (46 1 in dimensional form. and I p were transformed including the ref- 
erence values If, p s = p s c 2 s /^ and p s in the kernel and dummy variable of the 



integrals. Eq. (48) reads in dimensional notation 



rh l v l + p 



fh l e + v^p 



-{d-t+X) 







C.J 








2H 





PC. J - Po 



PC j - Po 



(48) 



It is interesting to note, that the system of equations representing the planar 



detonation ( 19 1 is, for a gas of known composition, a mono-parametric system 



dependent on D. The system represented by (48) depends not only on D, but 
also on curvature H and on stretch \. 

The velocity of the detonation D also suffers a change due to curvature that 
can be expressed through an asymptotic expansion of the form D = Dry + 
Dm)Z~ . The velocity of the detonation must be derived from the species 
equation. This derivation exceeds the scope of this paper and will be included 
in a forthcoming publication. 
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2.6. Virtual surface tension 



The normal momentum Eq. ( 46 1 is discontinuous across the jump with linear 
proportionality on curvature. Following the arguments in [11], an analogy can 
be established with an interface separating two immiscible fluids |15j in order 
to calculate the virtual surface tension of the detonation. At the interface, 

[p - <r nn ] =2Ha, (49) 

where a represents the surface tension. Our initial scale considerations suggest 
that the tangential stress is negligible. Allowing for mass transfer across the 
surface (e.g. due to evaporation) the formula must be modified 

[m 1 ^ 1 + p] =2Ha. (50) 



We may, identify terms with eq. ( 48 ) to obtain 



/ PCJ-P0f\ ~ ,p s c s pcj-Po , \ , r .. 

a = - J x - — Ip = -If [ I x —p s Ip = lfa . (51) 

1 Pcj-Po J V 7 Pc. j -Pa 1 



The coefficient of surface tension in an infinitely thin gas-dynamic discontinuity 
equivalent to a detonation is equal to the difference between the integral of the 
pressure between the CJ and the VN points minus the integral of the density 
between the same integration limits normalized by a factor. 

The evaluation of ao and the partial factors a i = —I^pgC^/'y and ao2 = 
(Pcj~Po)(pcj — /5o) _1 Ps^p is included in the Figure[3]for a gas of the indicated 
characteristics and different initial pressures po- The surface tension exhibits an 
inverse proportionality to the initial pressure. The existence of a minimum is 
due to the sum of the «oi an d a 02 factors (dashed lines) that combined creates 
the final dependency. 
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Figure 3. 



Figure [4] shows the dependence of the surface stress on the fuel concentration. 
An increase on the fuel concentration strongly increases the surface tension of 
the equivalent jump. 

Figure 4. 

The existence of virtual surface tension has strong implications for the sta- 



bility of the detonation. In this sense, Eq. (48) shows that the tangential 



momentum must not be continuous across the interface. The derivative of the 



surface tension, RHS of Eq. (48), has an analogous meaning to the Marangoni 
forces. A detonation may also propagate isotropically in a medium of constant 
composition forming cellular structures PQ, [12]. The RHS of the tangential 



impulse in Eq.(48) will cancel. The normal momentum component in Eq. (48) 
can be exploited to show that the mean curvature is constant at leading order 
for a complete cellular structure. However, the curvature will change in time 
through the derivative dependence of the energy jump condition. 



3. Conclusions 

In this paper we have derived a model which describes the detonation phe- 
nomenon as a gas-dynamic discontinuity. Our model was derived from the 
Euler equations. Modified Rankine-Hugoniot conditions link the fluid fields on 
the two sides of the detonation surface. The amplitudes of the jump depend 
of the precise location of the discontinuity related to the detonation structure. 
The position is defined by a density integral which yields vanishing of the excess 
surface mass, leading to the continuity of the mass flux across the flame. The 
continuity of the mass flux simplifies the conditions substantially and allows for 
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physical interpretation. An additional pressure jump proportional to the cur- 
vature appeared. This may be interpreted as virtual surface tension. The jump 
conditions obtained are valid for arbitrary chemistry. We plan to implement 
numerical codes for the computation of detonations based on our model. 
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Figure 1: Profiles of hydrodynamic (solid) and reactive (dashed) detonation models. Global 
view. Profile obtained with ZND theory coupled with rarefaction wave, for a gas of charac- 
teristics po = 100 kPa, po = lkg/m 3 , Q = 0.1 M J, 7 = 1.4, k = 1 10 5 s _1 , E/R g = 10000 K. 
designates normal status, VN von Neumann peak, CJ Chapman-Jouget point. 
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Figure 2: Profiles of hydrodynamic (solid) and reactive (dashed) detonation models. Detailed 
area. Profiles obtained with the same conditions as in Figure^ 
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Figure 3: Dependence of ag factor to the degree of over-driven detonation for a gas of char- 
acteristics po = lkg/m 3 , Q = 0.1MJ, 7 = lA,k = 110 5 s" 1 , E/R g = 10000 K obtained for 
different pressures 
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Figure 4: Dependence of ag factor to the degree of over-driven detonation for a gas of char- 
acteristics po = 100 kPa, p = 1 kg/m 3 , 7 = lA,k = 1 10 5 s _1 , E/R g = 10000 K obtained for 
different enthalpies of formation 
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Figure 1: Profiles of hydrodynamic (solid) and reactive (dashed) detonation 
models. Global view. Profile obtained with ZND theory coupled with rarefac- 
tion wave, for a gas of characteristics po = 1 00 kPa, po = 1 kg/m 3 , Q = 0.1 MJ, 
7 = 1.4, k = 1 10 5 s _1 , = 10000 A". designates normal status, ViV von 
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Figure 2: Profiles of hydrodynamic (solid) and reactive (dashed) detonation 
models. Detailed area. Profiles obtained with the same conditions as in Figure 

D 

Figure 3: Dependence of ao factor to the degree of overdriven detonation 
for a gas of characteristics po = 1 kg/m 3 , Q = 0.1 MJ, 7 = 1.4, k = 1 10 5 s _1 , 
E/R g = 10000 K obtained for different pressures. 

Figure 4: Dependence of «o factor to the degree of overdriven detonation 
for a gas of characteristics po = 100 kPa, po = 1 kg/m 3 , 7 = 1.4, k = 1 10 5 s _1 , 
E/Rg — 10000 K obtained for different enthalpies of formation. 
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